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ABSTRACT 

A fixed base helicopter simulation was implemented 
which may be easily altered to represent a variety of con- 
ventional single rotor aircraft. Test subjects performed 
90-second tracking tasks utilizing two cockpit displays for 
longitudinal control during landing approach. The performance 
data generated were compared with predicted values from an 
optimal pilot modeling technique. Tentative conclusions 
were drawn as to the effectiveness of each display as well 
as the predictive power of the pilot model. 



4 



TABLE OF CONTENTS 



I . INTRODUCTION 10 

II. METHOD OF INVESTIGATION 11 

III. APPARATUS AND PROCEDURES 11 

IV. OPTIMAL PILOT MODEL AND PREDICTIONS 15 

V. RESULTS 17 

VI. CONCLUSIONS 18 

APPENDIX A - THE DIGITAL PROGRAM 32 

APPENDIX B - THE ANALOG PROGRAM 40 

LIST OF REFERENCES . 42 

INITIAL DISTRIBUTION LIST 43 



5 



LIST OF TABLES 



I. UH-1H LONGITUDINAL DERIVATIVES AS USED IN 

SIMULATION 

II. TURBULENCE MODELS 

III. MODEL PERFORMANCE PREDICTIONS 

IV. SIMULATION PERFORMANCE DATA 



6 



LIST OF FIGURES 



1. Lines of Communication 24 

2. Apparatus 25 

3. Cockpit Display-- 26 

4. Optimal Pilot Model 27 

5. Performance Values, Subject A 28 

6. Performance Values, Subject B 29 

7. Control Motion, Baseline Display 30 

8. Control Motion, Quickened Display 31 



LIST OF SYMBOLS 



g 

I 

y 

M 

m 

q 

U o 



2 

Acceleration of gravity, ft/sec . 

Aircraft mass moment of inertia about y stability axis, 
slug-f t^. 

Moment about y stability axis, ft-lbs. 

Aircraft mass, slugs. 

Aircraft angular rate about y stability axis, radians/sec. 
Aircraft reference flight velocity, ft/sec. 



u Perturbation flight velocity along x stability axis, ft/sec. 

u Horizontal wind gust velocity, ft/sec. 
g 

w Perturbation flight velocity along z stability axis, 
f t/sec , 

Wg Vertical wind gust velocity, ft/sec. 

X Force component along x stability axis, lbs. 

Z Force component along z stability axis, lbs. 

6^ Cyclic pitch control input, ft. displacement at the 
pilot's hand. 

6 Collective pitch control input, ft, displacement at 
the pilot's hand. 

0 Perturbation pitch angle about y axis, radians. 



The following stability derivatives are defined for a 
reference flight condition in stability axes. 

Y = 1 3X 
q m 3q 

Y = 1 3X 
u m 3u 
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I . INTRODUCTION 



The use of vertical take-off and landing aircraft has 
increased significantly in both military and civil applica- 
tion. Helicopter flight under instrument flight rules has 
been accomplished on a routine basis for quite a while but 
low speed and terminal operations still require visual 
references outside the aircraft. The military services have 
had some success in low speed instrument flight but it is 
generally conceded that much work needs to be done in the 

I 

area of cockpit inforn ' One popular approach 



has been the use of electronic displays utilizing cathode- 
ray-tubes which have the advantage of almost unlimited 
flexibility in format and symbolism (References 3 and 2). 
Additionally, extensive work has been done in the field of 
pilot modeling and pilot-in-the-loop system performance (Ref. 3); 
but in order to optimize pilot performance it is necessary 
to consider the nature and extent of information presented 
to the pilot. This project was aimed at the construction 
of a fixed base helicopter flight simulator and the comparison 
of real pilot performance with pilot model predictions under 
a limited range of conditions and using two visual displays. 
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II , METHOD OF INVESTIGATION 



A helicopter simulation was constructed to represent 
the longitudinal flight dynamics of an unaugmented (i,e., 
without the usual mechanical stabilization bar) UH-1H 
helicopter. A fore and aft moving control stick provided 
cyclic pitch control inputs and a collective pitch level 
produced power commands. A hybrid computer gave real time 
solutions to the helicopter equations of motion, generated 
each of two cockpit displays, and computed performance 
data. Performance values for each display mode were 
analytically predicted using optimal pilot modeling tech- 
niques and compared with simulation data for two test 
subjects. 



The longitudinal motion of a helicopter can be described 
by the following equations of motion: 



III. APPARATUS AND PROCEDURES 





b c 



Z u - Z w 
u g w g 
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0 



q 



h = -w + U 0 • 

o 

6^ represents cyclic control input, 6 ^ collective input, 

and u and w horizontal and vertical wind gust velocities, 
g g 

Table I lists the stability derivative values used in this 
simulation. The equations above incorporate the following 
simplifying assumptions: 

1. The vehicle has been idealized as a rigid airframe to 
which is attached a rotor. 

2. The rotor is described by its tip path plane whose 
orientation determines propulsive and aerodynamic forces 
and moments. 

3. No rotor degree of freedom is considered other than 
control inputs which serve to describe instantaneous 
tip path plane orientation. 

4. Coupling between longitudinal and lateral motion has 
been ignored. 

5. Motion has been linearized as small perturbations about 
a predetermined reference flight path, 

The hybrid computer used consisted of a Scientific Data 
Systems SDS 9300 digital computer, A Comcor CI-5000 analog 
computer, and an Adage AGT-10 stroke written graphics 
processor, as illustrated in Figure 1. The digital computer 
was used to control the analog and graphics processor and 
to generate performance data, A program listing is found in 
Appendix A. In addition to the two display modes to be 
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described, provision has been made for the inclusion of 
flight director symbolism in a third display mode for future 
study. The analog computer gave solutions to the equations 
of motion in real time (Appendix B), and the graphics processor 
displayed textual and graphic information on a cathode ray 
tube in front of the test subject. Pilot control inputs 
were made directly to the analog using electrical voltage 
signals. The spring restrained cyclic stick operated a gear 
driven potentiometer and was linear in both force and dis- 
placement, The collective lever was unrestrained except for 
a friction brake and also used a gear driven potentiometer 
to generate a signal linear in displacement. The friction 
brake required a slight breakout force. Figure 2 shows the 
physical arrangement of the apparatus. 

Atmospheric turbulence was modeled by a sum of sine 
waves : 



5 

u = Z A.sinco.t, where 
g i=l 1 



and similarly, 

5 

w = Z A.sinw.t, where 
g i= l 1 i 



A ± = 4.472 

A 2 = 3.536 

A 3 = 2,236 

A 4 = 2,738 

A- = 2.236 
5 

A x = 4,472 
A 2 = 3.536 
A 3 = 2.236 
A 4 = 2.738 
A 5 = 2.236 



— 0 . 07 0 

w 2 = 0,210 

co 3 = 0.350 

a) . = 0 . 768 
4 

(l)r = 1.765 
o 

o)^ = 0.420 
u> 2 = 1-260 
a >3 = 2.100 
to 4 = 4 , 608 

u- = 10.050 

o 
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The amplitudes and frequencies above were chosen so that 
the integration over w of the power spectral densities 
closely approximated that for the turbulence spectra shown 
in Table II (Ref. 4). The advantage of this turbulence 
representation is that it is repeatable, has a mean value 
of zero, and appears completely random to the test subject. 

Two different displays were presented to the test 
subjects for evaluation. The first was basically a real 
world display similar to current instrument installations and 
consisted of a horizon bar and aircraft reference for pitch 
control along with a glideslope deviation indicator. Air- 
speed deviation was provided by angular displacement of the 
aircraft reference "wings." The second display had symbolism 
identical to the baseline display but was quickened with rate 
information according to the equations: 

0 . . = 0 + k *0 

quick 0 

h . , = h + k*h • 
quick h 

The quickening gains, k^ and k^, were simple pure gain 

approximations to the pilot transfer functions obtained in 
the modeling procedure. The display, shown in Figure 3, was 
6 1/2 inches wide and 7 inches high. The nominal eye-to- 
display distance was 2 1/2 feet. 

After considerable informal training, each subject 
received a formal training session in which he used the base- 
line and quickened display for 90 second runs each. In the 
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formal data sessions, each subject used each display for ten 
90 second runs. The data sessions went as follows: 

The baseline was flown for seven runs, followed by the 
quickened for seven. Then the quickened was flown for 
three runs followed by the baseline for three. Five of 
the "best" runs for each display were then selected. For 
ten data, runs with each display, mean scores were computed 
for each of the root mean square performance measures 
(u, 9,0, h, 6, , 6 ). Then each RMS score for each display 

l) C 

and run was normalized by dividing it with the respective 
mean RMS score. The normalized RMS scores were then summed 
to provide a single scalar index of performance for each 
display and run. 

The data in Table IV represent the means and standard 
deviations of the RMS scores for the runs with the five 
lowest indices of performance for each display. 

Figures 5 and 6 graphically portray the more pertinent 
performance data of Table IV for subject A (non-pilot) and 
subject B (pilot). Also indicated are the model predictions 
which best correlate with the data, i.e. F = 0,25. 

IV. OPTIMAL PILOT MODEL AND PREDICTIONS 

The analytic pilot model is shown schematically in 
Figure 4, and is basically the same as that offered by 
Kleinman and Baron (Ref. 5). This optimal pilot model 
includes the following assumptions: 
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1 . 



Subject to inherent limitations, the well trained and 
motivated pilot behaves in an optimal manner. 

2. With certain modifications, the solution to an optimal 
linear control and estimation problem will predict the 
pilot's control characteristics. These modifications 
are : 

a. A pure time delay is included in each of the pilot's 
control outputs. 

b. Each neuromuscular output system may be modeled as 
a first order lag. 

c. Each variable that the pilot observes is assumed to 
contain pilot induced additive white noise which 
scales with the variance of the observed variable. 

d. Each control output is assumed to contain pilot 
induced additive white noise which scales with the 
variance of the control. 

e. The pilot is able to perceive the first derivative 
of a variable, also noise contaminated, but no 
higher derivative. 

The index of performance for the optimization procedure was 
chosen subjectively by the analyst to represent the task 
objectives as perceived by the pilot. Model performance 
predictions are given in Table III for each of the displays. 

The parameter F represents the fraction of the pilot's 
attention devoted to the control task as a whole. For example, 
F = 0.25 suggests that 25% of the pilot's attention is devoted 
to control and 75% is dedicated to non-control activity such 
as communication, etc. 
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V . RESULTS 



Table IV lists the root mean, square performance data. 
Figures 5 and 6 give a graphical comparison of the most 
pertinent performance measures and the model predictions that 
best correlate with the data (F = 0.25). For each subject, 
the mean value is indicated, as well as plus and minus one 
standard deviation. In each case the pilot model predicted 
performance with a fair degree of accuracy. Additionally, 
the model predicted improvement for each performance measure 
using the quickened display, and the test subjects verified 
this prediction. However, the magnitude of performance 
improvement was greater in each case than was anticipated. 

It is possible that the task as perceived by the pilot may 
not be the same for each display mode. The task as understood 
may be to keep the deviation symbols within a finite distance, 
say 1/8", of their nominal position, and indeed there may 
exist a threshold of error below which the pilot makes no 
control input at all. Since the gains driving the deviation 
symbolism vary considerably between modes, each display then 
defines a new task. One possible arrangement would be to 
provide additional symbolism for the task variables, i.e. 
airspeed and glideslope deviations, which would be common 
between modes. 

The dramatic reduction in cyclic motion, 5 ^, was not 
predicted by the model. Typical time histories of cyclic 
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motion, Figures 7 and 8, show the nature of this improvement. 
Again the task defined by the display mode, in this case 
pitch attitude control, may be different. For both subjects, 
glideslope deviation control was poorest in comparison with 
model predictions, which suggests that more pilot attention 
was given to attitude/airspeed control than to the glideslope. 
Furthermore, collective motion was less for both subjects 
than was predicted. While the cyclic input was of a continuous 
nature, collective inputs were a series of steps (Figures 7 
and 8) in each case. Such control action probably represents 
a wait and see technique used by pilots as a result of the 
lag in collective response. This phenomenon may in part 
explain the discrepancy between predicted and measured 
glideslope deviations since the optimal control model did not 
address this type of control motion. 

Pilot comments indicated that the quickened display 
was more fatiguing than the baseline and represented an 
increase in pilot workload. Unfortunately, no measure of 
pilot workload, such as number of control inputs per minute, 
was recorded. Furthermore, no objective measure of the "well 
trained and motivated" assumption was available for each 
subject . 



VI , CONCLUSIONS 

The use of quickening in the cockpit display significantly 
reduced control motion and improved pilot performance. The 
predictive power of the optimal pilot model seemed adequate 
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even though this report represents only a first iteration. 
The quantitative quality of the model should improve as 
adjustments are made. The model index of performance can 
be adjusted to account for the pilot's apparent greater 
concern with pitch attitude/rate deviations as compared with 
altitude/rate deviations. In addition, the displays could 
be revised to eliminate any variation in the task presented 
to the pilot. As the control theoretic pilot model is 
refined and validated through testing it should become a 
valuable tool in the development of cockpit information 
systems and the expanded utility of the helicopter. 
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TABLE I 



UH-1H LONGITUDINAL DERIVATIVES AS USED IN SIMULATION 
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TABLE II 



TURBULENCE MODELS 



Vertical Turbulence 



$ (oj) 
w ' 

g g 



2a 2 L 
w w 

U 



L 

o 



ft 2 rad/sec^ 



o =5 ft/sec 
w ' 



L w = 100 ft 

U =60 knots (101.3 ft/sec) 
o ' 



Horizontal Turbulence 



Vg (u> 



2o 2 L 
u u 

“U 



L 0 
i+( 



ft 2 rad/sec 2 



a =5 ft/sec 
u ' 
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TABLE III 



MODEL PERFORMANCE PREDICTIONS 



F 


Baselin 

1.0 


e Display 
0.5 


0.25 


Qu 

1.0 


ickened I 
0.5 


)isplay 

0.25 


a u 

(ft/sec) 


.881 


1.50 


2,69 
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1.23 


2.45 


a 

w 

(ft/sec) 


2.94 
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CT © 

(rad/sec) 
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.0128 


.0231 


0 h 

(ft) 


1,81 


2,55 


4.01 


1,34 


2.01 


3.36 


a z 

6 B 

(ft) 
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.0300 


,0335 


.0256 


.0277 


.0311 


6 C 

(ft) 


,0513 
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.0675 


.0460 


.0511 


.0611 
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TABLE IV 






SIMULATION PERFORMANCE DATA 







Baseline 


Quickened 


Subject 


1 


2 


1 


2 


a 

u 


5 . 5St 


2.60 


3.91 


1.86 


(ft/sec) 


1 . 55tt 


.627 


.653 


.713 


a 0 


.0505 
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.0169 
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(rad/sec) 
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°0 
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a h 
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(ft) 
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Q 

o> 

O 
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.0292 


(ft) 
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.00712 
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.00501 



t mean 

ft standard deviation 
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Figure 1. 



Lines of communication. 





Figure 2. 



Apparatus . 
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Display Element 


Situation 


Units 


1. 


horizon bar 


pitch attitude 


20 degrees/inch 


2. 


airspeed indicaton 


airspeed error 


0.5 knots/degre 


3. 


aircraft symbol 


stationary 




4. 


glideslope bar 


glideslope deviation 


50 feet/inch 


5. 


glideslope indicators 


stationary 


+ or - 50 feet 



Figure 3. Cockpit display. 



26 





r 



~i 




z 

0 

I- 

< 

> 

01 
LjJ 
CO 
CD 
O 



UJ 

CO 



UJ 

O 

O 



O 

_) 

CL 





27 



Figure 4. Optimal pilot model. 
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Baseline display Q - Quickened display 
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Figure 6. Performance values, subject 
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Figure 7. Control motion, baseline display. 
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Figure 8. Control motion, quickened display. 



APPENDIX A 



DIGITAL PROGRAM 



INTEGER IGD(6),FRAME(20),GSLP(12),ACREF< 20) ,HORIZ(25) 
DIMENSION ITD ISO) , I TEXT! 12 ) ,UG( 1000) ,WG( 1000) 

NAMELIST MODE ,E16, El 7 ,E21 ,E23,E24,E25 
FRAME(l) = I HE AD ( 0 , 6 ) 

* 

DO 110 I =2 » 1 2 
READ ( 5 t 2 70) X,Y , I DM 
110 FRAME(I) = I PACK (X,Y, IDM) 

GSLP(l) = IHEAD( 0,7) 

« 

DO 120 1=2,11 
READ (5,270) X,Y, IDM 
120GSLPU) = IPACK( X ,Y , IDM) 

ACREF(l) = IHEAD( 0,9) 

* 

DO 130 1=2,19 
READ (5,270) X,Y,IDM 
130 ACREF(I) = IPACK (X ,Y, IDM) 

HORIZ(l) = I H EAD (0,8) 

-ajc 

DO 140 1=2,5 
READ (5,270) X, Y, IDM 
140 HORIZ(I) = IPACK(X,Y, IDM) 

A 1 =4. 472 
A 2=3. 536 
A3=2. 236 
A4=2. 733 
A5=2. 236 
CMU 1=0.07 
0MU2=0. 21 
0MU3=0 . 35 
0MU4=0 • 768 
0MU5= 1 . 765 
OMW 1=0. 42 
0 MW 2=1 . 26 
OMW 3= 2 . 1 
0MW4=4 . 608 
0MW5=10 .05 
UDR=9 .95 

£#:}::?:£:£ $:*:* if jjt £ ^ # # $####* If. # ## ### $ ♦ # # 

* TAKE CARE OF PERIPHERAL EQUIPMENT * 

OUTPUT ( 102) ' ENGAGE PATCHBOARDS , SELECT INPUT CONTROL, EX 
1ECUTE GATED, AND SET SENSE SWITCH* 

OUTPUT l 102) ' SELECT SETTINGS FOR POTS 016,017,021,023,0 
124,025* 

OUTPUT( 102) ' AS E 16 = ,E17= ,E21= ,E23= ,E24= ,E25= ' 
INPUT (101) 

* 

CALL SETPOT < 4H? 000 ,. 0091 , 4HP00 1 ,. 0053 , 4HP002 ,. 0044 ,4H 
1P004, .5680, 4HP00 5, .5395, 4HP006, .3000, 4HP007 ,.5067, 4HP0 

210. . 978 0, 4H P011 , . 1639, 4HP01 2, .0600,4HP0 13 , . 0 1 78 , 4 HPO 14 

3. . 4732.4HP015, . 1500,4HP020, .25 34,4HP022 , .6667 , 4HP027 , . 
47880, 4HP030, .642 0 , 4HP031 , . 3 5 52 , 4HP032 , . 6 800 , 4HP033 , . 15 
569, 4HP 037 , . 0 1 66 , 4HP 04 2 , . 3099 ,4HP045, . 1 4 0 0 , 4HP 05 0 , . 7725 
6 , 4HP0 5 1 , . 0242, 4HP052 , .0315, 4HP053 ,.l 173, 4HP054, . 7 077,4 
7HP056 , . 4000 , 4H P J 5 7 , . 3 500 , 4H P02 5 , E 25 , 4HP0 2 1 , E2 1 , 4H PO 23 , 
8E23,4HP024,E24,4HP01o,E16,4HP017,E17) 

IDEV = 1 

145 I F ( SE NS E S W I TC H ( 1 ) ) 155,150 

150 1F( SENSES WITCH( 2) ) 152,145 

152 I DE V= 2 
155 CONTINUE 



32 



CALL DGINIT { IDE V » I GD » 6 » l ER ) 

CALL DTINIT { IDEV , I T D » 6 0 , IER) 

160 OUTPUT ( 102) ' SELECT DISPLAY MODE, MODE=l OR 2' 

I NPUT ( 101 ) 

ENCODE (48,280, I TEXT) 

CALL TEXTO ( I DEV , I T EX T , 1 2 , 1 0 , 1 , 2 , 2 , I ER ) 

ENCODE (48,290, I TEXT) 

CALL TEXTO ( I DE V , I TEXT , 1 2 , 1 2 , 1 , 2 , 2 , 1 ER ) 

ENCODE (48 ,300, I TEXT ) 

CALL TEXTO ( I DEV , I TEXT , 12 , 14 , 1 1 2 , 2, I ER) 

ENCODE (48,310, I TEXT) 

CALL TEXTO ( I DE V , I TEX T , i 2 , 1 6 , 1 , 2 , 2 , 1 E R ) 

ENCODE (48,320, I TEXT) 

CALL TEXTO ( IDEV , I TEXT , 12 , 20 , 1 , 2 , 2 , 1 ER ) 

ENCODE (48,330, I TEXT) 

CALL TEXTO ( 1 DE V , I TE XT , 1 2 , 2 2 , 1 , 2 , 2 , 1 ER ) 

ENCODE (48,340, I TEXT) 

CALL TEXTO ( I DEV , I TEXT , 12 , 24, 1 , 2 , 2 , 1 ER) 

T = 0 • 0 

DO 163 1=1,1000 

SU1=0MU1*T 

SU2=OMU2*T 

SU3=0MU3*T 

SU4=0MU4*T 

SU5=0MU 5* T 

SW1=0MWI*T 

SW2 =0MW2*T 

SW3=0MW3*T 

S W4=0MW4*T 

S W5 =OMW 5* T 

UG( I) = ( 1./25.00) *( A1*SIN( SU1)+A2*SIN(SU2 )+A3*SIN< SU3) + 
1 A 4* SI N ( SU4)+A5*SIN( SJ5 ) ) 

WG ( I ) = ( l./25.00)*< A1*SIN( SW 1 ) + A2*S IN ( SW 2 ) + A3*S l N ( S W3 J + 
1 A4*SIN( SW4)+A5*SIN< SW5) ) 

T = T + ( 1 . O/UDR) 

163 CONTINUE 

165 IF (TEST(l).GT.O) GO TO 165 
CALL DT IN IT (IDEV, ITD, 60, IER) 

AB = 0. 0 

DO 170 1=6,10 
HORI Z ( I 1 = 0 
J = 1+7 
FRAME ( J ) = 0 
170 CONTINUE 

IF (MODE.GT.l) GO TO 200 



* MODE It BASIC DISPLAY 



£ 



CALL WRITECLOCK ( 0. ) 

CALL RESET (500) 

CALL COMPUTE 
CALL STARTCLOCK 

DO 180 1=1,7200 
CALL READCLOCK (V) 

IF(V.GT .5400.0) GO TO 190 
AB = AB+1 .0 
UG1 =UG( I ) 

WG1=WG( I I 

CALL DAC { 1 , UG1 , 2 , WG1 ) 

CALL ADK (0,AS, 1 , P , 2 , GS , 3 , Q , 4, GSDi 
P = . 2865* P 
GS = 0 . 4*G S 

Y = .1-P 

HORI Z ( 2 ) = I PACK ( - » 6 » Y » 0 ) 

HOP IZ ( 3 ) = IPACK (-0. 125, Y , 1 ) 

HOR 12(4) = IPACK< 0. 125, Y,0) 

HORI 1(5) = I PACK l . 6 , Y , 1 ) 

Y = -0.2-GS 
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GSLP<10) = I PACK (-0. 5 7 , Y , Q ) 
GSLP(ll) = IPACK (-0.37, Y, l ) 

GSLP( 12 ) = I PACK ( -0 . 5 7 , Y , 1 I 
ANGL = 0.7029* AS 
SAN = S IN ( ANGL ) 

CAN = COS( ANGL) 

Y = 0.1*SAN+0.1 
Y 1 = Y+0.015 
Y2 = Y-0.015 
Y3 = 0. 125*SAN+0.1 

X = -0. l*CAN-0. 125 

XI = X+0. 01 5*SAN 
X2 = X-C.015*SAN 

X3 = -0.12 5* CAN-0. 125 
ACREF! 2 ) = IPACK ( X » Y, 0 ) 

ACREFI3) = IPACK(X1,Y1,1) 

ACPEF(4) = I PACK (X3, Y3, 1 ) 

ACREF(5) = IPACK<X2,Y2. 1) 

AGREE! 6 ) = IPACK(X,Y,1) 

X = -X 

XI = -XI 
X 2 = -X 2 
X 3 = -X 3 

ACREFU5) = IPACK(X,Y,l) 

ACRE F( 1 6 ) = IPACK(X2,Y2, 1) 

ACRE F( 1 7 I = I PACK ( X3 * Y3 » 1 ) 

ACRE F( 1 8 ) = IPACK! XI, Y1 , 1 ) 

AC RE F( 19) = I PACK ( X , Y , l ) 

CALL GRAPHO ( I DEV , FRAME, 20, 1 , I ER ) 
CALL GRAPHO ( I DE V , GSL P , 1 2 , 2 , I ER ) 
CALL GRAPHO ( I DE V , ACk E F , 20 , 3 , I E R ) 
CALL GRAPHO ( I DE V , HOR I Z , 2 5 , 4 , I E R ) 
180 CONTINUE 
190 CALL HOLD (500) 

GO TO 260 



: # £ # 



'r 'r -v* “v 






* MODE 2, QUICKENED DISPLAY 

Vr Vr Vr ^ ^ ^ ^ ^ 'V -4# X X X X X X X X O* X X X X X X X V/ X X X X X X X X X X X 

V v -nr ^ -r t ^ ^ -nr* *r» nr ^ *r* -nr* -nr* X* -r -nr- nr nr ^ nr* nr nr n' n^ nr nr nr nr v ^r ^r 






"t* "T* V r V n' *Y' ^ 1* r 

'T 

j, o- ^ J> x ^ j, 

'(* *»- •v r A' t -r r *r t 



200 CALL WRITECL0CK (0.) 

CALL RESET (500) 

CALL COMPUTE 
CALL STARTCL0CK 

DO 210 1=1,7200 
CALL READCLOCK ( V) 

IF! V.GT. 5400.0) GO TO 220 
AB = A8 + l • 0 
UG1 =UG( I ) 

WG1 =WG ( I \ 

CALL DAC ( 1 , UG 1 » 2 , WG 1 ) 

CALL ADK (0,AS, 1 , P , 2 , GS » 3 , Q , 4 , G SD ) 
P = . 2665*! P+2. 760*. 4*Q) 

GS = 0 . 4* ( GS + . 2*GSD ) 

Y = . 1-P 

HOR I Z ( 2 ) = I P ACK ( - • 6 , Y , 0 ) 

HOR I Z ( 3 ) = IPACK(-0.125,Y, 1 ) 

HORI Z ! 4 ) = IPACKtO. 125,Y,0) 

HOR IZ ( 5 ) = IPACK! .6, Y, 1) 

Y = -0.2-GS 

GSLP(IO) = IPACK(-0.57,Y,0) 
GSLP(ll) = IPACK(-0.37,Y, 1) 

GSLP! 12 ) = IPACK(-0.57,Y, 1) 

ANGL = 0 . 7029*A S 
SAN = SIN! ANGL) 

CAN = COS ( ANGL ) 

Y = 0. l*SAN+0. 1 
Y 1 = Y+0.015 

Y 2 = Y-0.015 

Y3 = 0. 12 5*SAN + 0 . 1 

X = -0. l*CAN-0. 125 

XI = X+0.015*SAN 
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X 2 = X-0.015*SAN 
X 3 = -0. 125*CAN-0. 125 
ACREF ( 2 ) = I P ACK ( X , Y » 0 ) 

ACREF( 3 ) = I PACK t X 1 , Y 1 , 1 ) 

ACR E F( 4 ) = IPACK(X3,Y3, 1) 

ACREF( 5 ) = IPACK (X2,Y2, 1 ) 

ACREF(6 ) = I PACK ( X , Y , 1 ) 

X = -X 

XI = -XI 
X2 = -X 2 
X3 = -X 3 

ACRE F( 1 5 ) = IPACK(XtYil) 

ACREF( 16) = IPACK(X2, Y2, 1) 
ACREF(17) = I P AC K ( X3 , Y 3 , 1 ) 

ACREF (18) = I PACK ( X 1 , Y 1 » 1 ) 

ACPEF( 19) = IPACK(X, Y, 1) 

CALL GRAPHO ( l DE V , FRA ME , 20 , 1 , I E R ) 
CALL GRAPHC ( I DE V , GSL P , 1 2 , 2 , I E R ) 
CALL GRAPHO ( I DEV , ACREF ,20, 3, IER) 
CALL GR A D HO ( IDE V , HOR I Z , 2 5 , 4 , I E R ) 
210 CONTINUE 
220 CALL HOLD (500) 

260 CONTINUE 



* THIS SECTION TAKES INTEGRATED SQUARE VALUES FROM THE 

* ANALOG AND GIVES ROOT MEAN SQUARE PERFORMANCE VALUES. 

3^ Jfc jfc ^ ix jJs jJc iis Jjt 3jx r r ^ ^ ^ A y 



CALL DG IN IT ( I DE V , IGD , 6 , I ER ) 

CALL ADK( 16,SAS, 1 7 , SG S , 20 , S DC , 2 1 , SDB , 22 , THD , 23 , THE ) 
CALL RESET (500) 



SAS = 


SAS/ ( .104* 


E25 ) 


SGS = 


SGS/( .012* 


E21 ) 


SDC = 


SDC / ( 480 . * 


E23 ) 


SDB = 


SD9/ (120.* 


E24 ) 


THD = 


THD / ( 360 . 0 


* E 1 7 ) 


THE = 


THE/( 2250. 


0*E 16 ) 


SAS = 


ABS( SAS ) 




SGS = 


A BS ( SGS) 




SDC = 


ABS(SDC) 




SDB = 


ABS( SDB) 




THD = 


A BS ( THD ) 




THE = 


ABS ( THE ) 




RMAS 


= SQRT(SAS) 




RMGS 


= SQRT(SGS) 




RMDC 


= SQRT(SDC) 




RMDB 


= SQRT(SDB) 




RMTD 


= SQRT(THD) 




RMTH 


= SQRT { THE ) 




UDR = 


AB/120.0 





WRITE (6,350) MODE 
WRITE (6,360) UDR 

WRITE (6,370) RMA S , RMG S , RMDB ,RMDC , RMTD , R MTH 
GO TO 160 

270 FORMAT (2F1 
280 FORMAT (• 

1 ') 

290 FORMAT (• 

1 * ) 

300 FORMAT (• 

1 • ) 

310 FORMAT (' 

1 • ) 

320 FORMAT ( • 

1 ') 

330 FORMAT (' 

1 ' ) 

340 FORMAT ( * 

1 • ) 

350 FORMAT («1» 



0.4,11) 

THIS IS A HELICOPTER TRACKING 
PROBLEM REQUIRING THE ADJUSTMENT OF 
POWER AND PITCH ATTITUDE TO MAINTAIN 
60 KTS ON A 6 DEGREE GLIDESLOPE. 

WHEN READY TO BEGIN, PRESS THE 
RED BUTTON ON THE COLLECTIVE LEVER. 
THE TASK WILL LAST FOR 90 SECONDS. 

* 'DISPLAY MODE ',11/) 
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H- * -) 



360 FORMAT {'O't'UPDATE RATE AVERAGED ',F8.5,' TIMES PER S 
1EC0ND' / ) 

370 FORMAT < 'O', 'RMS AIRSPEED ',F8.5,» FT/SE 

1 C ' / ' 0 ' » 'RMS HEIGHT ',F8.5, FT'/'O', 

2 'RMS CYCLIC *,F8.5,' FT'/'O', 

3 'RMS COLLECTIVE ',F8.5»' FT'/'0' f 

4 'RMS VERTICAL T UR8ULENCE *tF8.5,' FT/SEC'/'O', 

5 'RMS PITCH RATE *,F8.5,' RAD/ SEC ' / ' 0 ' , 

6 'RMS PITCH ANGLE ',F8.5,‘ RAD '/) 

END 
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PROGRAM VARIABLES 



AB 

ACREF 

ANGL 

AS 

CAN 

E16 

E17 

E21 

E23 

E24 

E25 

FRAME 

GS 

GSD 

GSLP 

HORIZ 

IDEV 

IDM 

IER 

IGD 

ITD 

I TEXT 

MODE 

P 



Incremental number of display updates during simulation. 
Vector of coordinates defining aircraft reference. 

Angle defining airspeed deviation. 

Airspeed deviation. 

Cosine of given angle. 

Setting for potentiometer 016. 

Setting for potentiometer 017. 

Setting for potentiometer 021. 

Setting for potentiometer 023. 

Setting for potentiometer 024. 

Setting for potentiometer 025. 

Vector of coordinates defining display frame. 

Glideslope (height) deviation. 

Glideslope deviation rate. 

Vector of coordinates defining glideslope reference. 
Vector of coordinates defining horizon reference. 
Graphics terminal in use, 1 or 2. 

Graphics instruction, draw or move. 

Error type code. 

Graphics instruction, size of graphics input. 

Graphics instruction, size of text input. 

Vector of text input. 

Display mode, 1 or 2. 

Pitch angle. 
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Q 


Pitch rate. 


RMAS 


Root mean square airspeed deviation. 


RMDB 


Root mean square cyclic movement. 


RMDC 


Root mean square collective movement. 


RMGS 


Root mean square glideslope deviation. 


RMTD 


Root mean square pitch rate. 


RMTE 


Root mean square pitch angle. 


SAN 


Sine of given angle. 


SAS 


Mean square airspeed deviation. 


SDB 


Mean square cyclic movement. 


SDC 


Mean square collective movement. 


SGS 


Mean square glideslope deviation. 


T 


Time increment. 


THD 


Mean square pitch rate. 


THE 


Mean square pitch angle. 


UDR 


Update rate for the 90 second period. 


UG 


Vector of horizontal turbulence magnitudes. 


UG1 


Specific horizontal turbulence value. 


WG 


Vector of vertical turbulence magnitudes. 


WG1 


Specific vertical turbulence value. 


V 


Internal clock count, 60 cycles per second. 


X 


Running variable, abscissa. 


Y 


Running variable, ordinate-. 


A1,A2, . 


..A5 Sine wave magnitudes. 



0MU1 ,OMU2 , . . , OMU5 Horizontal turbulence frequencies 
0MW1 , 0MW2 , , . , 0MW5 Vertical turbulence frequencies. 
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SU1 , SU2 , , 
SW1,SW2, . 



. , SU5 
. , SW5 



Sine values for horizontal turbulence. 
Sine values for vertical turbulence. 
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APPENDIX B - ANALOG PROGRAM 
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